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ELECTROMAGNETIC FIELD ANALYZER, ELECTROMAGNETIC FIELD 
ANALYZING PROGRAM, AND STORAGE MEDIUM FOR RECORDING THE 

PROGRAM 

TECHNIC A fcBACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates to electromagnetic field 
simulation technology used in development developing and 
designing high-frequency products such as multilayer 
products, giga filters, and EMI (Electro-Magnetic 

Interference) elimination filters, etc . , — and more particularly^ 
the present invention relates to an electromagnetic field 
analyzer, an electromagnetic field analyzing program, and a 
storage medium for recording the program in which the 
simulation time can be ghortcnod is reduced . 
Background 

2. Description of the Related Art 

In recent years, the development and designif^ of high- 
frequency products have been vigorously performed aggressively 
pursued, and, in order to perform an electromagnetic field 
analysis of high-frequency products, electromagnetic field 
simulation softwares have been widcly sof tware has been 
commonly used. Generally, there are many cascs ins t ances where 
a finite element method is used for electromagnetic field 
analysis. In the electromagnetic field simulation using the 
finite element method, since it takes a long time to solve 
simultaneous linear equations, high-speed solving methods are 



desired. In particular, in the field of electromagnetic field 
analysis, since nothing except — for only a direct method called 
a Gauss 1 elimination method can be used, the problem of talcing 
a— long calculation time times is more serious. 

However, by^a paper (R. Hiptmair, "Multigrid mcthod Method 
for Maxwell's cquat iono Equat ions , " SIAM Journal of Numerical 
Analysis, vol. 36, no. 1, pp. 204-225, 1 999-)- ) , which is 
referred to herein as Non-patent document 1 , and which was 
made public by Hiptmair in 1999, it is proved that a solution 
method called a multi^grid method can be also used in an 
electromagnetic field analysis and the method io paid 
attention to as a greatly speeding up solution when reduces 
the calculation time as compared with to a direct method. 

In the multi^grid method, an analysis object is divided 
into two classes of fine elements and coarse elements^ and the 
solution of fine elements is found by making use of determined 
using the solution of coarse elements. 

When simultaneous linear equations are calculated by using 
the Gauss' elimination method, generally the calculation time 
is proportional to the third power of the dimensions of a 
matrix. Since the size of the dimensions corresponds to the 
number of elements, when the number of elements becomes 
doublc is doubled , the calculation time becomes is 8 times as 
large, and, when the number of elements becomes is 1 0 times as 
large, the calculation time becomes is 1^000 times as large. 
Then, — when When the Gauss ' direct method is applied to only 
coarse elements, an approximate solution to fine elements is 
obtained by using a— the solution to coarse elements, and the 
Gauss' elimination method is not used, the calculation time 

2 



for fine elements can be reduced to 1/8 and 1/1000, 
respectively. Since additional calculation time is added, 
although the practical calculation time is not reduced a-s— to 
such a great extent , the calculation time can bc is always 
greatly reduced without fail . 

The multi^grid method is classif icd can be broken down into 
two methods: a method using a nested mesh and a method using 
a non-nested mesh. However, since a prolongation matrix is 
incomplete as described later, there is nothing except — §ef 
using a the method using a nested mesh must be used . 

Figs. 1 2A and 1 2B arc drawings for describing show the 
division into elements using a nested mesh. Fig. 1 2A shows 
the division into elements using a nested mesh in the case of 
two-dimensional analysis where a coarse triangular element is 
uniformly divided to form fine triangular elements. 
Furthermore, Fig. 1 2B shows the division into elements using a 
nested mesh in the case of three-dimensional analysis where a 
coarse tetrahedral element is uniformly divided to form eight 
fine tetrahedral elements. As -ts— understood from Figs. 1 2A 
and 12B, there are geometrical restrictions between coarse 
elements and fine elements in the case of when using a nested 
mesh. 

Nor— Another prior art method is described in non- patent 
document 1 : R . — Hiptmair , — - 2: D. Dibben and T. Yamada, "Non- 
Nested Multigrid method for Maxwell's equations, " SI AM Journal 
of Numerical Analysis, — vol . — 3-6-? — Ne^ — — pp~. — 204 225, — 1 900 

Non patent document 2: — B-z — Dibben and T. — Yamada, — "Non 
nested multigrid and automatic mesh coarsening for high 
frequency electromagnetic problcmo and Automatic Mesh 
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Coarsening for High Frequency Electromagnetic Problems " IEEJ 
Investigating Research Committee matcrial Material , SA-02-34, 
pp. 71-75, 2002^ 

Figs. 1 3A and 1 3B arc drawings — for dcscribing de scribe a 
first problem e# that occurs when using a nested mesh. Fig. 
1 3A shows how a circle is divided into elements using a nested 
mesh to form coarse elements. Furthermore, Fig. 1 3B shows how 
coarse elements are divided into fine elements using a nested 
mesh. As irS— under stood from Figs. 1 3A and 13B, a curved 
surface cannot be correctly expressed because of the 
geometrical restrictions of a nested mesh. 

Furthermore, Fig. 1 3C shows how a circle is divided into 
elements using a non-nested mesh to form fine elements. As irs- 
understood from Fig. 13C, since there are no geometrical 
restrictions when using a non-nested mesh , a curved surface 
can be correctly expressed. 

Figs. 1 4A to 1 4C arc to describe a second problem o# that 
occurs when using a nested mesh. Fig. 1 4A shows how a square 
is divided into elements by— using a nested mesh to form coarse 
elements. Furthermore, Fig. 1 4B shows how coarse elements are 
divided into fine elements. As irs— understood from Figs. 1 4A 
and 14B, coarse elements cannot be partially divided into fine 
elements because of the geometrical restrictions of a nested 
mesh ." 

Furthermore, Fig. 14C shows how a square is divided into 
elements by— using a non-nested mesh to form fine elements. As 
-i-s— understood from Fig. 14C, coarse elements can be partially 
divided into fine elements because e# there are no geometrical 
restrictions when using a non-nested mesh . 
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Most industrial products have curved surfaces , such as irft 
a circular cylinder and a— square surfaces, and they inevitably 
contain partially fine portions. Accordingly, it is difficult 
to apply a nested mesh to industrial products. On the other 
hand, since there is are no geometrical restrictions in a non- 
nested mesh, fine elements can be freely formed. However, 
since the positional rclation relationship is not systematic, 
it is difficult to make the electromagnetic field of coarse 
elements related to the electromagnetic field of fine elements 
It is a prolongation matrix to make coarse elements related to 
fine elements, but no precise prolongation matrix has been 
found. This means that a multi^grid method cannot be 
practically used, and accordingly, nothing can be used except 
for a Gauss' direct method^ which takes much requires a lot of 
time^ in the analysis of magnetic field. 
Disclosure of Invention 

The present invention has been done to solve the 

SUMMARY OF THE INVENTION 

To overcome the problems described above problem, — and it 

is an object , preferred embodiments of the present invention 
^e— provide an electromagnetic field analyzer, an 
electromagnetic field analyzing program, and a recording 
medium having the program recorded t here in t her eon , in which 
the analysis of an electromagnetic field can be made by 
using per formed via a multi^grid method using a non-nested mesh 

It is another object of the present invention to provide ^ 

and an electromagnetic field analyzer, an electromagnetic 
field analyzing program, and a recording medium in which the 
program has been recorded, in which the time required for an 



electromagnetic field analysis is greatly reduced. 

According to an aspcct a preferred embodiment of the 
present invention, an electromagnetic field analyzer 
compriocG includes dividing means for dividing form data as an 
analysis object into coarse elements and fine elementst^ 
forming means for forming a matrix making de fining an 
electromagnetic field vector of the coarse elements divided by 
the dividing means related to an electromagnetic field vector 
of the fine elements, and calculating means for calculating an 
approximate solution of the electromagnetic field vector of 
the fine elements by applying an iteration method of 
simultaneous linear equations while referring to the matrix 
formed by the forming means. 

Preferably, the forming means forms the matrix by 
expressing the elements of electromagnetic field vectors at 
the sides at a fine element using an interpolation function in. 
the coarse elements. 

More preferably, when the length of a side li of the fine 
element is |ljj, an interpolation function showing the 
rclation relationship between the electromagnetic field at 
position x of the side li of the fine element and the 
electromagnetic field at a side j of the coarse element is Nj C 
(x) , and the tangential vector of the side li of the fine 
element is ti, the forming means (25) — forms a matrix Pij 
preferably by Expression (14) to be described later. 

More preferably, in the calculation calculating means, the 
accuracy of an approximate solution of the electromagnetic 
field vector of the fine elements is improved such that high- 
frequency components included in the approximate solution of 
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the electromagnetic field vector of the fine element are 
removed by applying a stationary iteration method of 
simultaneous linear equations, such that a residual in the 
fine elements is mapped to a residual in the coarse elements 
by using the matrix formed by the forming means, such that a 
correction vector to the coarse elements is formed by applying 
a direct method or a non-stationary iteration method of 
simultaneous linear equations, and such that a correction 
vector to the fine elements is obtained by using the matrix 
formed by the forming means. 

According to another aopcct pref erred embodiment of the 
present invention, in a computer program to cauoc control a 
computer to execute an electromagnetic field analyzing method 
for analyzing an electromagnetic field as an analysis object, 
the computer contains includes a first storage means for 
otoring device arranged to store divided elements and a second 
storage means for otoring device arranged to store a matrix, 
and the electromagnetic field analyzing method 
contains includes the steps of dividing form data as an 
analysis object into coarse elements and fine elements and 
storing the elements in the first storage mcans device , forming 
a matrix making def ining an electromagnetic field vector of the 
coarse elements stored in the first storage mcQno device 
related to an electromagnetic field vector of the fine 
elements and storing the matrix in the second storage 
mcanG device , and calculating an approximate solution of the 
electromagnetic field vector of the fine elements by applying 
an iteration method of simultaneous linear equations while 
referring to the matrix stored in the second storage 
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mcans device . 

According to another aspcct pref erred embodiment of the 
present invention, in a computer readable storage medium— 
being computer readable , — for recording having stored thereon a 
program to causc control a computer to execute an 
electromagnetic field analyzing method for analyzing an 
electromagnetic field as an analyzing object, the computer 
contains includes a first storage means for storing device 
arranged to store divided elements and a second storage means 
for storing device arranged to store a matrix, and the 
electromagnetic field analyzing method contains includes the 
steps of dividing form data as an analysis object into coarse 
elements and fine elements and storing the coarse and fine 
elements in the first storage mcans device , forming a matrix 
making de fining an electromagnetic field vector of the coarse 
elements stored in the first storage mcans device related to an 
electromagnetic field vector of the fine elements and storing 
the matrix in the second storage mcans device , and calculating 
an approximate solution of the electromagnetic field vector of 
the fine elements by applying an iteration method of 
simultaneous linear equations while referring to the matrix 
stored in the second storage mcans device . 

According to an aspect of the present invention, — since 

Since the forming means forms a matrix malcing de fining an 
electromagnetic field vector of coarse elements divided by the 
dividing means related to an electromagnetic field vector of 
fine elements, it became is possible to perform an 
electromagnetic field analysis by using via a multi^grid method 
using a non-nested mesh. Furthermore, since the calculating 
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means calculates an approximate solution of the 
electromagnetic field vector of fine elements while referring 
to the matrix formed by the forming means, ar^ the calculating 
means is not required to calculate the electromagnetic field 
of fine elements by using a direct method-? — it became possible 
to largely reduce a . Thus, the time necded required for an 
electromagnetic field analysis is greatly reduced . 

Furthermore, since the forming means forms a matrix such 
that an electromagnetic field vector element at a side of fine 
elements is expressed by using an interpolation function in 
coarse elements, it became possible to a matrix is easily form 
a matrix f ormed . 

Furthermore, in the forming means, since a matrix Pij is 
formed preferably by Expression (14) to be described later, a 
matrix can bc is even more easily formed. 

Furthermore, in the calculating means, an approximate 
solution of an electromagnetic field vector of fine elements 
is corrected by using a matrix formed by the forming means-? — i-fe- 
became possible to improve . Thus, the accuracy of the 
approximate solution of an electromagnetic field vector of 
fine elements is greatly improved at high speed. 
Brief Description of the Drawings 

Other features, elements, steps, characteristics and 

advantages of the present invention will become more apparent 
from the following detailed description of preferred 
embodiments thereof with reference to the attached drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is an external appearance of shows an 
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electromagnetic field analyzer according to a ^a preferred 
embodiment of the present invention. 

Fig. 2 is a block diagram showing the construction of the 
electromagnetic field analyzer according to aft a preferred 
embodiment of the present invention. 

Fig. 3 is a block diagram showing the functional structure 
of the electromagnetic field analyzer. 

Fig. 4 is a flow chart for describing the procedure of the 
electromagnetic field analyzer. 

Fig. 5A shows a microstrip line divided into large 
elements . 

Fig. 5B shows a microstrip line which is analyzed using 
the large elements shown in Fig. SAj_ and in which the elements 
are subdivided based on their error estimation. 

Fig. 6 shows components of an electromagnetic field 
which that are in contact with the sides of a tetrahedron. 

Fig. 7 shows the rclation relat ionship between an 
electromagnetic field vector u F in Q F and an electromagnetic 
field vector u c in fi c . 

Fig. 8 shows one example of a rectangular waveguide as an 
analysis object of an electromagnetic field analyzer according 
to -etPt a preferred embodiment of the present invention. 

Fig. 9 shows the calculation time to the number of 
elements when the rectangular waveguide shown in Fig. 8 is 
driven in TE1 0 mode of 2.45 GHz, which is compared between the 
case where a related direct method is used and the case where 
a method of a preferred embodiment of the present invention is 
used . 
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Fig. 10 shows one example of a patch antenna as an 
analysis object of the electromagnetic field analyzer 
according to etna, preferred embodiment of the present invention. 

Fig. 11 shows the calculation time to the number of 
elements when the patch antenna shown in Fig. 10 is driven at 
7.0 GHz, which is compared between the case where a related 
direct method is used and the case where a method of the 
present invention is used. 

Fig. 1 2A shows elements divided by using a nested mesh in 
the case of two-dimensional analysis. 

Fig. 1 2B shows elements divided by using a nested mesh in 
the case of three-dimensional analysis. 

Fig. 1 3A shows coarse elements formed such that a circle 
is divided into elements by using a nested mesh. 

Fig. 1 3B shows fine elements into which coarse elements 
are divided by using a nested mesh. 

Fig. 1 3C shows fine elements formed such that a circle is 
divided into elements by using a non-nested mesh. 

Fig. 1 4A shows coarse elements formed such that a square 
is divided into elements using a nested mesh. 

Fig. 1 4B shows fine elements into which coarse elements 
are divided by using a nested mesh. 

Fig. 14C shows fine elements formed such that a square is 
divided into elements by using a non-nested mesh. 
Reference Numerals 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Fig. 1 computer main body 
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2 display device 

3 FD drive 

4 FD 

§ keyboard 

6 mouse 

7 CD ROM device 

8 CD ROM 

9 network communications device 

W GP¥ 

W ROM 

V± RAM 

Vi hard disk 

— — — ~2A form data storage portion 

-2-2 division into element portion 

2^ clement storage portion 

£4 vector/ matrix — forming portion 

2-5 prolongation matrix forming portion 

2-6 vector/matrix storage portion 

2-7 prolongation matrix storage portion 

— approximate value calculation portion 

2-9 approximate solution correction portion 

3-6 electromagnetic field vector storage portion 

Best Mode for Carrying Out the Invention 

Fig- — ^ — is an external appearance of shows an 

electromagnetic field analyzer according to a-na preferred 
embodiment of the present invention. The electromagnetic 
field analyzer contains pref erably includes a computer main 
body 1 , a display device 2, an FD (Flexible Disk) drive 3 to 
which an FD 4 is attached, a keyboard 5, a mouse 6, a CD-ROM 
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(Compact Disk- Read Only Memory) drive 7 to which a CD-ROM 8 
is attached, and a network communications device 9. An 
electromagnetic field analyzing program is supplied through a 
recording medium such as an FD 4, or a CD-ROM 8, etc, or other 
suitable media. An electromagnetic field analysis is 
performed such that the electromagnetic field analyzing 
program is executed by the computer main body 1 . Furthermore, 
the electromagnetic field analyzing program may be supplied to 
the computer main body 1 by ant her another computer through the 
network communications device 9. 

Fig. 2 is a bock diagram showing the construction of the 
electromagnetic field analyzer of an according to a preferred 
embodiment of the present invention. The computer main body 1 
shown in Fig. A — contains 1 includes a CPU (Central Processing 
Unit) 10, a ROM (Read Only Memory) 11, a RAM (Random Access 
Memory) 12, and a hard disk 13. The CPU 10 conduct s per forms 
processing while the CPU 10 inputs and outputs data to the 
display device 2, the FD drive 3, the keyboard 5, the mouse 6, 
the CD-ROM device 7, the network communications device 9, the 
ROM 11, the RAM 12, or the hard disk 13. The electromagnetic 
field analyzing program in the FD 4 or the CD-ROM 8 is 
preferably stored i^on the hard disk 1 3 through the FD drive 3 
or the CD-ROM drive 7 . The electromagnetic field analysis is 
performed such that the CPU 1 0 properly loads the 
electromagnetic field analyzing program from the hard disk 13 
in the RAM 12 and executes the program. 

Fig. 3 is a block diagram showing the functional structure 
of the electromagnetic field analyzer according to -afta 
preferred embodiment of the present invention. The 
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electromagnetic field analyzer con tains includes a form data 
storage portion 21 for storing form data as an analysis object, 
a division-into-element portion 22 for dividing the form data_^_ 
such that the size of elements is changed, an element storage 
portion 23 for storing the elements divided by the division- 
into-element portion 22, a vector /matrix forming portion 24 
for forming a vector and matrix in simultaneous linear 
equations, a prolongation matrix forming portion 25 for 
forming a prolongation matrix, a vector/matrix storage portion 
26 for storing the vector and matrix formed by the 
vector /matrix forming portion 24, a prolongation matrix 
storage portion 27 for storing the prolongation matrix formed 
by the prolongation matrix forming portion 25, an approximate 
value calculation portion 28 for calculating an approximate 
value of an electromagnetic field vector, an approximate 
solution correction portion 29 for correcting the approximate 
solution calculated by the approximate value calculation 
portion 28, and an electromagnetic-field vector storage 
portion 30 for storing^ as an electromagnetic field vector^ 
the approximate solution after having been corrected by the 
approximate solution correction portion 29. 

Moreover, the form data storage portion 21, the element 
storage portion 23, the vector /matrix storage portion 26, the 
prolongation matrix storage portion 27, and the 
electromagnetic field vector storage portion 30 are provided 
in a fixed area in the RAM 12 or on the hard disk 13 shown in 
Fig. 2. Furthermore, the function f unctions of the division- 
into-element portion 22, the vector /mat rix forming portion 24, 
the prolongation matrix forming portion 25, the approximate 
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value calculation portion 28, and the approximate solution 
correction portion 29 is rcalizcd are arranged such that the 
CPU 10 shown in Fig. 2 executes the electromagnetic field 
analyzing program loaded into the RAM 12. 

Fig. 4 is a flow chart for describing the procedure of the 
electromagnetic field analyzer according to e^ a preferred 
embodiment of the present invention. First, the division- 
into-element portion 22 reads form data as an analysis object 
stored in the form data storage portion 21 , divides the form 
data into elements, and pre serves stores the divided elements 
in the element storage portion 23 (S11). The division-into- 
element portion 22 divides the form data into triangular or 
quadrangular elements in the case of two-dimensional analysis, 
and divides the form data into tetrahedral, triangle-pole- 
shaped, or square-pole-shaped elements in the case of three- 
dimensional analysis . 

The division-into-element portion 22 divides the form data 
while changing the size of the elements. In order to make the 
description simple, here elements of two different sizes are 
formed and it is assumed that a set of larger (coarse) 
elements is made fi c and a set of smaller (fine) elements is 
made Q F . Moreover, a superscript of c represents coarse, and a 
superscription of F represents fine. 

The sets of Q c and Q F may be independently formed, 
bu^t however , when the set of £2° is first formed and then the set 
of Q F is formed by subdividing the elements, the processing 
time is shortened. reduced. Regarding the method of 
subdivision, there arc is a method for subdividing elements 
which are larger in volume and a method for analyzing elements 
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by using Q c and subdividing elements having _a__large error. 
Either of them may be used, b^ however , the method in which 
error is evaluated has a— better calculation accuracy. 
Moreover, regarding the method for evaluating error, Reference 
Document No. 1 (No'rio Kamiya et al . , "computer and 
application ' 3 monk, — special edition/error estimation in 
software and adaptive clement," No . "Computer and Application's 
Monk, Special Edition/Error Estimation in Software and 
Adaptive Element," No. 42, CORONA PUBLISHING CO., 1993 )) can j_ 
may be referred to. 

Figs. 5A and 5B show examples of division into elements in 
a microstrip line. In Fig. 5A, the microstrip line is divided 
into larger elements. Furthermore, in Fig. 5B, an analysis is 
performed by using the larger elements in Fig. 5A and the 
elements are subdivided based on their error estimation. 
Since the electromagnetic field is concentrated at the edge of 
an electrode, the portion is subdivided into smaller elements. 

Next, the vector/matrix forming portion 24 reads the sets 
£1° and Q F stored in the element storage portion 23, 
discretization of the sets takes place using their elements 
such that a finite element method is applied to a Maxwell 
equation, and a matrix and column vector in simultaneous 
linear equations are formed (S12). Here, the Maxwell equation 
as an object is represented by the following expression-s : 
Expression 1 

Vxff = (d + ja>i)E -0) 

VxE = -ja)fiH -(2) 

Here, Expressions 1 and 2 

Vx// = (ct + jO)i)E -(1) 
VxE = -jo)mH -(2) 
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Where E: electric field, H: magnetic field, e: dielectric 

constant, \i: permeability, o: conductivity, co: angular 

frequency, V: differential operator (d/dx, d/dy, d/dz) , X: 

outer product, j: imaginary number 

From these expressions, the following differential 

equation can bc is obtained— j_ 

Expression -3-_3 
Vx pVxu-O) 2 qu-0 ■■■(3) 

u represents an electric field, magnetic field, etc., and 
p and q represent characteristic values of dielectric constant 
and permeability. For example, when Expression (3) is 
obtained by substituting Expression (2) into Expression (1) 
and eliminating H, u represents an electric field, p represent 
the reciprocal of a permeability, and q represents a 
dielectric constant . 

When a finite element method is applied to an element V k 
formed in the above former process described above , the 
following expression can be is obtained—: 

Expression 3 

1 {V*N u )p{Vxu)dV -a> 2 j N ki qudV = [ y {N u x P Vxu)ndS 

-(4) 

Here , — represents the surface of and n represents an 

outward normal unit vector. represents — an interpolation 

function — (also called a shape — function, — insertion — function , 

basis function, — etc . ) When tctrahcdral elements and linear 

equations arc used, — the interpolation function is given by the 
following expression . 
Expression 4 
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l k (V x N ki )p (V x u) dV - 1» 2 {, A^fa/P = (N» x £V x «> ndS 

-(4) 

Here, 3V k represents the surface of V k and n represents an 

outward normal unit vector. N k j represents an interpolation 
function (also called a shape function, insertion function, or 
basis function) . When tetrahedral elements and linear 
equations are used, the interpolation function is given by the 
following expression : 



Expression 5 



^1 = 








\1 2 \{W7L,- 


-Avl 3 ) 


N k3 = 


|/ 3 |(AVi 2 - 




* k4 = 






N kS = 






N k6 = 







Here, L x to L 4 represent a volume coordinate system and 
I 111 to I lei represent the length of a side lj.. Moreover, 
regarding the detail of the volume coordinate system, 
Reference Document 2 (T. Itoh, G. Pelosi and P. P. Silvester, 
"Finite Element Software for Microwave Engineering," John 
Wiley & Sons, pp. 101-125, 1996) easmay be referred to. 

Fig. 6 shows components of an electromagnetic field u 
which that are in contact with the sides of a tetrahedron. The 
components of the electromagnetic field u are disposed 
arranged so as to be disposed along the six sides, and they 
are represented by ui to U6- The electromagnetic field at 
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point x inside the element V k can bc is represented by the 

following expression using Expression (-5-6) and the six 

components of an electromagnetic field— 

Expression 5 

6 

i=i 

When Expression — (-6-) — is substituted into Expression — f4-) — etf*d 

applied to all the elements of Q ? - f — the following simultaneous 
linear equations can be obtained. 
Expression 6 

«(*) = £ N u {x)u u -(6) 

Here , — a f is vectors which arc electromagnetic fields 

disposed at the sides, — and the dimensions arc equal to the 

number of the sides as a general rule, The components of a 

matrix A ¥ are formed by the following expression. 

When Expression (6) is substituted into Expression (4) and 

applied to all of the elements of Q F , the following 
simultaneous linear equations are obtained: 
Expression 7 
A F u F = b F -(7) 

Here, u F represents vectors which are electromagnetic 

fields disposed at the sides, and the dimensions are equal to 
the number of the sides as a general rule. The components of 
a matrix A F are formed by the following expression: 
Expression 8 

III [L <?xN u )p<yxN v )dV-a> 2 l N u qN v dV] -(8) 

Furthermore, a vector b F is formed by the following 
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expression . 

Expression &9_ 

III LWm*??***)'"* "(9) 
Here, N F is the number of elements in f2 F . 

The vector/matrix forming portion 24 also forms a matrix 
A c corresponding to Q c -±** via the same procedure and preserves 
A F , b F , and A c in the vector/matrix storage portion 26. When u F 
is f ound determined by calculating simultaneous linear 
equations in Expression (7) , the electromagnetic field 
distribution can be calculated using Expression (6) to result 
iftproduce the solutions of Maxwell equations in Expressions 
(1 ) and (2) . 

Next, the prolongation matrix forming portion 25 reads 
elements Q F and Q c formed in the former process from the 
element storage portion 23 and, based on the geometric 
information, the prolongation matrix forming portion 25 forms 
a prolongation matrix P which malccs relates the electromagnetic 
field vector u F in Q F related to the electromagnetic field 
vector u c in Q c (S13). When the prolongation matrix is used, 
as shown in the following expressions^ one vector can make be 
used to determine the other vector determined . j _ 
Expression 9 



Expressions 1 0 and 1 1 
u F =Pu c -(10) 
u c =P'u F -(11) 

Here, t means a transposed matrix. 

Fig. 7 shows the relat ion relationship between an 
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electromagnetic field vector u F in Cl F and an electromagnetic 
field vector u c in £2 C . The j component Uj C of the column vector 
u c is disposed at the side of an element of £l c and the 
electromagnetic field u(x) at the position x in Q F can be 
represented using an interpolation function N c in Q c — _^ 
Expression 1 0 



j 

The i component U j. p of u p is made an average value of an 



electromagnetic — field at — the side 1 ± where — t&e — i component u ± 

is disposed . That is , — the i component u ^ F of u F is represented 

by the following expression. 
Expression 1 1 



Here , — fe-j. represents a tangential unit vector of the side 

i*-: When Expression — (12) — is substituted for Expression — (13) 

and compared with Expression — (10), — the i j — component ~P ±j of a 
prolongation matrix P is given by the following expression. 
Expression 1 2 

"(*) = £ N*j{x)u c j -(12) 

j 

The i component Uj F of u F is an average value of an 

electromagnetic field at the side lj where the i component Uj - 
is disposed. That is, the i component U j F of u F is represented 
by the following expression: 
Expression 1 3 

«f = lfrt«(*K<tf -(13) 

Here, tj represents a tangential unit vector of the side 
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lj. When Expression (12) is substituted for Expression (13) 
and compared with Expression (10), the ij component P j j of a 
prolongation matrix P is given by the following expression. 

Expression 14 

P v =^[Nj(xndl -(14) 

The prolongation matrix forming portion 25 forms a 
prolongation matrix using Expression (14) and preserves the 
prolongation matrix in the prolongation matrix storage portion 
27. 

Next, the approximate value calculation portion 28 sets an 
approximate initial value v F to u F in order to calculate the 
electromagnetic field vector v F using simultaneous linear 
equations according to a multi^grid method (S14). The initial 
value v F may be zero. Moreover, the det ail details of the 
calculation of the simultaneous linear equations according to 
the multij^grid method ea ftmay be referred to in Reference 
Document 3 (W. Briggs, V. Henson and S. McCormick, "A 
Multigrid Tutorial," SIAM) . 

Next, the approximate value calculation portion 28 removes 
high-frequency components of error contained in v F by applying 
a stationary iteration method of simultaneous linear equations^ 
such as the Jacob method, the Gauss-Seidel method— or the SOR 
method, etc . , to Expression (7) (S15). Moreover, the 
dctail details of the iteration method can be referred to in 
Reference Document 4 (Satomi Hasegawa, Hidehiko Hasegawa, 
Kiyoji Fujino, "Iteration Method Templates", Asakura 
shotcn Shoten , 1996). 

Next, the approximate value calculation portion 28 maps a 

residual r F in ft F to a residual r c in Q c by using the following 
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cxprcosion expressions in order to remove low- frequency 



components (S16). 



Expression 1 3 




-(15) 
-(16) 



Next, the approximate value calculation portion 28 
f inds det ermines solutions of simultaneous linear equations in 
the following expression (17) and calculates a correction 
vector e c (SM^U_ 

Expression 441 7 



In Expression (17), although a non-stationary iteration 
method^ such as a direct method according to the Gauss 
elimination method, a CG method — or a GMRES method-? — etc . , is 
preferably used, since Expression (17) is a small matrix 
whcn as compared with to Expression (7) , the calculation can be 
performed in a short time. Moreover, the dctail details of the 
iteration method can be referred to in the above -mentioned 
document 4 . 

Next, the approximate solution correction portion 29 
improves the accuracy of the approximate solution by returning 
the correction vector of Cl c to the correction vector of Q F by 
using the following Expressions (18) and (19) at (SI 8-hr) : 

Expression 4-& 1 8 and 1 9 




-(17) 



e F =Pe c 

F F 



F 



•••(18) 
•••(19) 
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The approximate solution correction portion 29 makes high- 
frequency components attenuated by applying smoothing shown in 
step 15 a couple of times to the approximate solution (S19) . 

Although the above processing in step S1 5 to S19 is the 
basic algorithm of a multi^grid method, since there are cases 
where the convergence is relatively slow, the convergence is 
quickened by accelerated by using another iteration method for 
simultaneous linear equations together as required (S20) . 
Concretely Specif ically , the above -described processing is used 
as a pretreatment for a non- stationary iteration method^ such 
as a CG method, a GMRES method-? — or a GCR method, etc . , — or the 
processing is accelerated by a residual cutting method, etc. 
Moreover, the dctail det ails of the residual cutting method can 
be referred to documcnt in Reference Document 5 (Kikuchi et al . , 
JSME International Journal, 62-604, Series B p. 4076—, 
December 1 996-4-3-). 

Next, the approximate solution correction portion 29 
evaluates the residual of Expression (15) in step (S21 ) . When 
the residual is large (S21 is no) , since the electromagnetic 
field vector does not converge, the process goes back returns 
to step S1 5 and the following processing is repeated. 
Furthermore, when the residual i^s- has a sufficiently small 
value (S21 is yes) , it is undcrstood determined that the 
electromagnetic field vector converges, and the approximate 
solution v F is stored in the electromagnetic field vector 
storage portion 30 (S22) and the processing is finished. 

Fig. 8 shows one example of a rectangular waveguide as .an 
analysis object of an electromagnetic field analyzer according 
to eft a preferred embodiment of the present invention. The 
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dimensions of this rectangular waveguide ieare as shown in Fig 
8 and the inside is filled with air. The electromagnetic 
field analysis is performed using tetrahedral primary elements 

In Fig. 9, the calculation time to the number of elements 
when the rectangular waveguide shown in Fig. 8 is driven in 
TE10 mode of 2.45 GHz is compared between the case where a 
related direct method is used and the case where a method e# 
according to the present invention is used. As irs— understood 
from Fig. 9, when the number of elements increases, the effect 
of the case where an electromagnetic field analyzer 
e #according to the present preferred embodiment is used is 
eminent more pronounced and, when the number of elements is 
206,266, the speed of the processing is about 17 times as high 
as when compared with greater than that of the related direct 
method. 

Fig. 10 shows one example of a patch antenna as an 
analysis object of the electromagnetic field analyzer 
according to a^ a preferred embodiment of the present invention 
The dimensions of the patch antenna are a^— shown in Fig. 10^ 
and the patch antenna is f ormed disposed on a dielectric 
substrate having a thickness of about 0.794 mm and a 
dielectric constant of about 2.2, The electromagnetic field 
analysis is performed by using primary tetrahedral elements. 

In Fig. 11, the calculation time to the number of elements 
when the patch antenna shown in Fig. 10 is driven at 7.0 GHz 
is compared between the case where a related direct method is 
used and the case where a method according to a preferred 
embodiment of the present invention is used. As is understood 
from Fig. 11, when the number of elements increases, the 



effect of the case where an electromagnetic field analyzer of 
the present preferred embodiment is used is cmincnt pronounced 
and, when the number of elements is 88,445, the speed of the 
processing is about 1 5 times as high as when compared with the 
than that of the related direct method. 

As described above, according to the electromagnetic field 
analyzer of the present preferred embodiment , since a 
prolongation matrix P making de fining the electromagnetic field 
vector u F in Q F related to the electromagnetic field vector u c 
in Q° is formed and an electromagnetic field vector is 
calculated by using this prolongation matrix, it became is 
possible to perform an electromagnetic field analysis by using 
a multi^_grid method using a non-nested mcch mesh . 

Furthermore, since it became is possible to perform an 
electromagnetic field analysis by using a multi^grid method 
using a non-nested mesh, the time required for analysis can be 
is greatly reduced when as compared with to an electromagnetic 
field analysis using a direct method. 

The embodiments described this time arc exemplification in 

all aspects and not restrictive. The scope of the present 

invention is not shown by the above description, but shown by 
the scope of the claims, — and the range equivalent to the scope 
of the claims and all modifications within the scope of the 
claims arc to be included. 

— While the present invention has been described with respect 
to preferred embodiments, it will be apparent to those skilled 

in the art that the disclosed invention may be modified in 
numerous ways and may assume many embodiments other than those 
specifically set out and described above. Accordingly, it is 
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intended by the appended claims to cover all modifications of 
the present invention that fall within the true spirit and 

scope of the invention. 



27 



